Dzyaloshinskii-Moriya Interaction and Spiral Order in Spin-orbit Coupled Optical 

Lattices 



Ming Gong 1 , Yinyin Qian 1 , V.W. Scarola 2 , and Chuanwei Zhang 1 '^3 

1 Department of Physics and Astronomy, Washington State University, Pullman, WA, 99164 USA 
2 Department of Physics, Virginia Tech, Blacksburg, Virginia 24061 USA 
3 Department of Physics, the University of Texas at Dallas, Richardson, Texas, 75080 USA 

We show that the recent experimental realization of spin-orbit coupling in ultracold atomic gases 
can be used to study different types of spiral order and resulting multiferroic effects. Spin-orbit 
coupling in optical lattices can give rise to the Dzyaloshinskii-Moriya (DM) spin interaction which 
is essential for spin spiral order. We derive an effective spin model in the deep Mott insulator region 
at half filling, and demonstrate that the DM interaction in optical lattices can be made extremely 
strong with realistic experimental parameters. The rich phase diagrams of the effective spin model 
for fermion and bosons are obtained via classical Monte Carlo simulations. 
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The interplay between ferroelectric and ferromagnetic 
order in complex multiferroic materials presents a set of 
compelling fundamental condensed matter physics prob- 
lems with potential multifunctional device applications 
Ferroelectric and ferromagnetic order compete and 
normally cannot exist simultaneously in conventional ma- 
terials. Coexistence of these orders often relies on strong 
correlation. In some strongly correlated systems, such 
as the perovskite transition metal oxides [2|-l7j , these two 
phenomena can occur simultaneously. Construction and 
design of high-T c magnetic ferroelectrics is still an open 
and active area of research [8[ . These and other complex 
solids incorporate different types of interactions, includ- 
ing electron-electron interactions, electron-phonon inter- 
actions, spin-orbit (SO) couplings, lattice defects, and 
disorder, making the determination of multiferroic mech- 
anisms a remarkable challenge for most materials. Due 
to the inherent complexity of the different interactions 
involved, the underlying mechanism for the emergence of 
multiferroic behavior is not fully clear 0, [To[ . An unbi- 
ased and direct method to explore multiferroic behavior 
(e.g., dynamics, critical exponents, the impact of con- 
trolled disorder, etc.) in an ideal setting is thus highly 
appealing. 

The realization of a superfluid to Mott insulator transi- 
tion of ultracold atoms in optical lattices [n[ opened fas- 
cinating prospects [12| for the emulation of a large vari- 
ety of novel magnetic states and other strongly correlated 
phases found in solids. In optical lattices different types 
of interactions can in principle be independently tuned 
and controlled, thus providing an ideal platform for the 
emulation of quantum phases. This platform can offer 
great insight into the physics of complex solids that incor- 
porate strong correlation. As a result, possible magnetic 
phase diagrams have been explored in the optical lattice 
setting [Til- 15]. For instance, it has been shown 13|, 1^ 
that the effective Hamiltonian of spin- 1/2 atoms in op- 
tical lattices is the XXZ Heisenberg model in the deep 
Mott insulator regime. On the experimental side, su- 



perexchange interactions between two neighboring sites 
have already been demonstrated [16| and quantum simu- 
lation of frustrated classical magnetism in triangle optical 
lattices has also been realized [17j. These experimental 
achievements mark the first steps towards the quantum 
simulation of possible magnetic phase transitions in op- 
tical lattices. 

In this letter, we show that the power of optical lat- 
tice systems to emulate magnetism can be combined with 
recent experimental developments [18l-l2l( realizing SO 
coupling to emulate multiferroic behavior. Recently, SO 
coupled optical lattices have been shown to possess fiat 
energy band 22], and important experimental progress 
[23| has been made for observing the single particle band 
dispersion (including the flat band) in SO coupled optical 
lattices. The main findings of this work are the follow- 
ing: (I) We derive the effective Hamiltonian for spin- 1/2 
fermions and bosons in optical lattices in the large in- 
teraction limit. We show that SO coupling leads to an 
effective in-plane Dzyaloshinskii-Moriya (DM) term, an 
essential ingredient in models of spiral order and multi- 
ferroic effects in general. We estimate that the DM term 
is of the same order as the Heisenberg coupling constant, 
J. (II) We study the phase diagram of the effective spin 
model using classical Monte Carlo (MC). We find that 
competing types of spiral order depend strongly on the 
SO coupling strength and effective Zeeman field. (Ill) 
We find that the critical temperature for the spiral order 
can be made rather high, ~ J. Thus, if magnetic quan- 
tum phase transitions can be emulated in optical lattices, 
then spiral order and multiferroics can also be realized in 
the same setup with the inclusion of SO coupling. 
Effective Hamiltonian: We consider spin-1/2 ultracold 
atoms loaded into a two-dimensional (2D) square optical 
lattice. We restrict ourselves to the deep Mott insula- 
tor regime where the charge/mass degree of freedom is 
frozen while the spin degree of freedom remains active. 
Here the atomic hyperfine levels map onto effective spin 
states. The scattering length between the atoms in op- 
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FIG. 1: Transition processes due to spin-conserving tunnel- 
ing (solid lines) and SO coupling mediated tunneling (dashed 
lines) for spin-1/2 fermions (a) and spin-1/2 bosons (b). The 
lowest 4 levels correspond to nearly degenerate ground states. 
The higher energy levels are the excited states. 



tical lattices can be controlled by a Feshbach resonance 
[24} . To derive the inter-spin interaction in this regime 
we first consider a two-site tight-binding model, 

H = - ^ t a c\ a C 2o . + Vso + i ^ U ™' : n i<T n i<r' : > (!) 

where c\ a creates a particle (either a boson or a fermion) 
in a Wannier state, wi )CTl localized at a site i and in 
a spin state a e {t,|}. n ia = c\ a c i(T and U a(T > = 
9aa' J dx.\wi, a \ 2 \wi :a >\ 2 , where g GG > is the interaction 
strength between species a and a' . Here :: denotes nor- 
mal ordering. For a general theory the tunneling is as- 
sumed to be spin dependent, which is a feature unique to 
ultracold atom systems [3, EH . The second term is the 
SO coupling, which may take various forms depending on 
the configuration of the laser fields [25]. Here we use the 
Rashba SO coupling for simplicity, 



V so = i\c\e z • (a x d)c • + h.c. 



(2) 



where c\ = (cL, cj^), d = (dx, dy) is the vector from site j 
to z, and A is the SO coupling strength. The magnitude 
and sign of A can be tuned in experiments [26[ using 
coherent destructive tunneling methods. 

In the deep Mott insulator regime, the effective Hamil- 
tonian can be obtained using the standard Schrieffer- 
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Let 



Wolf transformation 
V = — [iS, Hq] to eliminate the first-order term, we have 



the second order Hamiltonian H e ^ = Hq + \[iS, V]. In 
the spin representation we define Si = ^ ss c\ s (j SS 'C isf1 
and extend the two-site model to the whole lattice, yield- 
ing 

H ^ = E E JaSfs; + E B - s * + 

^D V -(S, \S, i • s, -r,, -s,, (3) 

The first two terms denote Heisenberg exchange and Zee- 
man terms, respectively, while the last two terms arise 
from SO coupling. In solid state systems the third term 
is called the DM interaction [28], which is believed to 
drive multiferroic behavior. In the following we derive 
the coefficients in the effective Hamiltonian by consider- 
ing the coupling between four internal degenerate ground 
states \a) e {| t;t), I t;l), I l;t), I Ui)} through the spin 
independent and dependent tunnelings t a and A. The 
couplings are different for fermions and bosons, as illus- 
trated in Fig. [TJ 

Fermionic atoms: For fermionic atoms, there are only 
two possible excited states |ex) = | ti? 0) an d |0;ti) 5 as 
shown schematically in Fig. [1] (a). The only relevant 
two-body interaction is = U. The couplings between 
the ground states \a) and the excited state |ex) are de- 
noted by solid lines for pure spin-conserving tunneling 
(t a terms) and dashed lines for SO coupling mediated 
tunneling (A terms). We find (J x + J y )/2 = U^t^/U, 
(J x - J y )/2 = S(-dx 2 + dy 2 )UX 2 /(U 2 - (/i t - /i;) 2 ), 
and J z = 2(t t - ttf/U - 4d 2 UX 2 /(U 2 - h 2 ), with 
d 2 = dx 2 + dy 2 . The DM interaction coefficient is 
D = 2(t t + h)(h 2 - 2U 2 )\/(U(h 2 - U 2 ))(dy,dx,0), and 
the effective Zeeman field B = 4((/i^- ji^)-2d 2 h\ 2 j (h 2 - 
/7 2 ))(0,0, 1). Note that without SO coupling the model 
reduce to the well-known XXZ Heisenberg model with 
rotation symmetry [3, fl4j . However, this symmetry is 
broken by the SO coupling, yielding an XYZ-type Heisen- 
berg model. Similar results are also observed for bosons. 
Bosonic atoms: For bosonic atoms, there are six ex- 
cited states |ex) = | tt;0), | t|;0), | ||;0),|0;tt>, |0;tl), 
|0;||), as shown in Fig. U (b). Without SO coupling, 
the only allowed inter-state second-order transition is be- 
tween 1 2) and |3), similar to the fermionic case. The pres- 
ence of SO coupling permits other inter-state transitions, 
therefore the bosonic case is much more complex than the 
fermionic case. For simplicity we only show the results 
for Uf\- = Un = U^i = U, which yields (J x + J y )/2 — 
-tt^tjU, (J x - J y )/2 = 4(d 2 x - d 2 y )UX 2 /(U 2 - h 2 ), J z = 
-4t t t l /U + 2((h 2 - U 2 )(t t -t l ) 2 + 2U 2 d 2 \ 2 )/U(U 2 - h 2 ), 
D = -2(t t + t^ih 2 - 2U 2 )X/U(h 2 - U 2 )(dx, dy, 0), and 
B = (0,0,4(/i t -/^)). 

The last term in Eq. [3] reads as • • Sj = 

V ^T00y (SfSy + SySf), where rj = +1(-1) for fermions 
(bosons). This term arises from the coupling between 
states |1) and |4), |1)(4| = SfSj -S» S» +i(SfS» +S» Sj). 
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Here the real part contributes asymmetric terms to the 
Heisenberg model, while the imaginary part contributes 
to Yij. In a square lattice with d x d y = 0, this term van- 
ishes. However, for a tilted lattice, such as triangular and 
honeycomb lattices, this term should be significant. 
Lattice parameters: We estimate the possible parameters 
that can be achieved in a square optical lattice V(x, y) = 
V(x)V(y), where V(x) = Vq sm 2 (kLx) = sEr sm 2 (kLx). 
The recoil energy Er = h 2 k 2 j /2m, where m is the mass 
of the atom, and fc^ is the wavevector of the laser. The 
SO coupling in free space is V so = ^(px^y —Py&x), where 
7 ~ hkR/m, kR is the wavevector of the external Raman 
lasers, and kR ~ fc^ in most cases. The Raman lasers 
are pure plane waves, and serve as a perturbation to the 
hopping between adjacent sites. 

We use the Wannier functions of the lowest band 
without SO coupling to calculate the tight binding pa- 
rameters: t = J dxw* (x)[p 2 1 '2m + V(x)]wi+i(x) and 
A = J dxw*(x)V so Wi+i(x). In a square lattice, coor- 
dinates decouple and the Bloch functions are Mathieu 
functions. The Wannier functions can be obtained from 
the Fourier transformation of the Bloch functions. Our 
numerical results are presented in Fig. [2] (a). The large s 
limit, t ~ t appr . = 4£^^/ v /7T8 3 / 4 exp(— is also plot- 
ted for comparison. Note that U/Er ~ (8/7r) 1 / 2 A:a s< s 3 / 4 
is in general much larger than t and can be controlled 
through Feshbach resonance independently. 

In Fig. [2] (b) we plot |D|/| J| as a function of f = X/t 
for U a(T > = U, t(j = t. |D|/|J| reaches the maximum value 
of 1.0 at A = t. This is in sharp contrast to models of 
weak multiferroic effects in solids with D / J ~ 0.001 — 0.1, 
which is generally induced by small atomic displacements 
[29j| . Optical lattices, by contrast, can be tuned to exhibit 
either weak or strong DM terms. This enhanced tunabil- 
ity enables optical lattice systems to single out the effects 
of strong DM interactions and study the impact of the 
DM term. 

There are notable differences between our model and 
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FIG. 2: (Color online), (a) Tunneling amplitudes as a func- 
tion of lattice depth, t is the hopping due to the kinetic 
energy, £ apP r. is the analytical expression in the deep lattice 
regime, and A is the SO mediated hopping strength, (b) Plot 
of |D|/|J| as a function of X/t for U aa > — JJ^ta— t. 
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FIG. 3: Phase diagrams of 2D optical lattices with fermions 
(a) and bosons (b) obtained from classical Monte Carlo sim- 
ulations for an 8 x 8 lattice at temperature T = 0.05 Jo- The 
phases diagrams are determined by the magnetization order, 
the spiral order, and the spin structure factor. Different re- 
gions correspond to: M = 0, P — for green, P — 
for grey, M = 0, P / for cyan, and M / 0, P / for red. 
The abbreviations are: (a) AF: antiferromagnetic phase with 
zero total magnetization; MAF: antiferromagnetic phase with 
non-zero total magnetization; NS: zero magnetization spiral 
order; MS: magnetic spiral order; NMFS: nonmagnetic flux 
spiral phase; MFS: magnetic flux spiral phase. In (b), SM: 
simply magnetic order; SMS: simply magnetic spiral order: 
Other abbreviations are the same as in (a). The dashed lines 
are guides to the eye. The spin structure factors of the points 
marked by plus signs are shown in Fig. [H 



corresponding models of solids (i) In solids the SO cou- 
pling arises from intrinsic (atomic) SO coupling and D 
is generally along the z direction (out of plane). How- 
ever, in our model D is in the plane and the out of plane 
component is zero, (ii) In our effective spin model, J°- 
depends on the direction of the bond (d X: d y ) and the SO 
coupling strength, while in solids J-*- is independent of 
SO coupling due to its negligible roles. 
Spiral Order and Multiferroics in 2D Optical Lattices: 
We explore the rich phase diagrams of the effective spin 
Hamiltonian using classical MC simulations. The classi- 
cal MC has been widely used to explore the phase dia- 
grams of the Heisenber g m odel with DM interactions in 
the context of solids [8|, l30l432j (and hence weak DM in- 
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FIG. 4: Spin structure factors for different quantum phases 
marked by plus sign in Fig. [3] The upper panels show the 
results for fermion at h/Jo = 1.1, while the lower panels show 
the results for boson at h/Jo = 0.218. 
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FIG. 5: (a) The spin configuration of fermions in an 8 x 8 
lattice at T — 0.05 Jo, £ = 1.0 and h/Jo = 1.5. The corre- 
sponding magnetization and spiral order as a function of tem- 
perature is shown in (b), with the inset plot Xv ~ (SP) 2 /T vs 
temperature, which indicates a phase transition at T c ~ 0.5 Jo. 
Similar features can also been found for bosons with the same 
parameters. 



teractions). This method may not be used to determine 
the precise boundaries between different phases but can 
be an efficient tool to determine different possible phases. 
Due to the unique features of our effective model (e.g., 
strong DM interactions) the phase diagrams we present 
here are much more rich and comprehensive than those 
explored in the context of solids. We focus on the regime 
where t a = t and U aa > = U (spin independent), and 
define Jo = 4t 2 /U as the energy scale. The rescaled ef- 
fective Hamiltonian becomes 

ff = ^ E j'SfSf + n-SixSj + h^s?, (4) 

ij a=x,y,z i 

where j* = -1 + (<% - ^)^ 2 , j y = -1 - (<% - ^)^ 2 , 
jz = -1 + £ 2 , D = -2£{dy, d x , 0), and £ = A/ Jo. 

Eq.[4]hosts a variety of magnetic and spin spiral phases, 
which are generally characterized by the magnetic and 
spiral order parameters |33i l34j 

m = J2 s i and p = d ^ x Si x ( 5 ) 

i (iJ) 

However, these two order parameters do not fully char- 
acterize the phase diagrams because in some cases there 
are still local magnetic or spiral orders although both M 
and |P| are vanishingly small. In these cases, we also 
take into account the spin structure factor: 

5(k) = ^(Si • Sj) exp(zk • (R* - R,)), (6) 

ho 

which shows peaks at different position in momentum 
space for different phases. For instance, the peak of 
the spin structure factor is at k = (0,0) for ferromag- 
netic phases, k = (7r,7r) for antiferromagnetic phases, 
and (7r,0) (or (0,7r)) for the flux spiral phase (P = 
but with nontrivial local spin structure). General spiral 
orders correspond to other k. We obtain the phase dia- 
grams by analyzing both the order parameters and spin 



structure factors. We have not checked for long range 
order in the spin structure factor. We expect quasi- long 
range order to accompany magnetized phases at low ft, 
e.g., a ferromagnetic phase for £ <C 1. 

Fig. [3] shows the phase diagrams of an 8 x 8 lattice 
for both fermions and bosons. The results are obtained 
after 10 6 thermalization steps followed by 10 6 sampling 
steps in each MC run at low temperature (T = 0.05 Jo). 
We have checked that for lower temperatures the phase 
diagrams do not change quantitatively. We also verify 
that similar phase diagrams can be obtained for larger 
system sizes, however, the spiral orders in larger opti- 
cal lattice become more complicated, and the boundary 
between different quantum phases is shifted. 

The phase diagrams show a rich interplay between 
magnetic orders and spin spiral orders. For instance, 
for fermions with small SO coupling (£ < 0.25), the 
ground states are anti-ferromagnetic states with zero 
(nonzero) magnetization for a Zeeman field h/Jo < 0.8 
(h/Jo > 0.8). While for large SO coupling (f > 1.45), 
the ground states are either nonmagnetic or magnetic 
flux spiral phases (similar to the flux phase with a small 
spiral order P). For £ ^> 1 the DM term is not important 
because D/J ~ l/£ ; therefore the pure flux phase with 
zero spiral order can be observed. Similarly, the increas- 
ing SO coupling for bosonic atoms gives rise to a series of 
transitions from simply magnetic (ferromagnetic at small 
h) order to simply magnetic spiral order (with zero total 
spiral order but local spiral structure), then to magnetic 
spiral orders (or non-magnetic spiral orders) and finally 
to flux spiral orders. The emergence of the spiral order 
and flux order with increasing SOC can be clearly seen 
from the change of the spin structure factors in Fig. [4j 
which shift from k = (0,0) or (7r,7r) to (-zr, 0) and (0,7r). 

The spin spiral order phase transition temperature is 
comparable to the magnetic phase transition tempera- 
ture, ~ Jo- In Fig. [5] (a), we plot the spin configura- 
tion of fermions at T = 0.05 Jo, £ = 1.0 and ft = 1.5 
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(MS phase), which show clear spiral ordering. The cor- 
responding order parameters P and M are plotted in 
Fig. [5] (b) as a function of temperature. The inset show 
susceptibility Xp ~ (SP) 2 /T. We see a phase transition 
at T c ~ 0.5 Jo, which is comparable to the magnetic crit- 
ical temperature [14| (in 2D, the Heisenberg model has a 
critical temperature T c = Jo in mean- field theory). Note 
that spiral order can also exist in the frustrated model 
without SO coupling, however, the critical temperature 
is generally much smaller than the magnetic phase tran- 
sition temperature [H, [35[. Our results therefore show 
that SO coupling in the absence of frustration provides 
an excellent platform for searching for high-T c spiral or- 
ders and multiferroics in optical lattices. 

Finally we note that different spiral orders may be ob- 
served using the optical Bragg scattering methods (36| . 
which probe different spin structure factors for different 
spiral orders. Similar methods have been widely used 
in solid state systems. Furthermore, in optical lattices, 
the local spin magnetization at each lattice site (thus the 
magnetic order M) as well as the local spin-spin correla- 
tions (thus the spiral order P) can be measured directly 
(37l . , which provides a powerful new tool for under- 
standing the physics of spiral orders and multiferroic ef- 
fects in optical lattices. 

Note Added. During the preparation of this manuscript 
we became aware of preprints 39|-41] on similar topics. 

Acknowledgement MG thanks S. Liang for numerical 
assistance with classical MC simulations. This work is 
supported by AFOSR (FA9550-11-1-0313), DARPA-YFA 
(N66001-10-1-4025, N66001-11-1-4122), ARO (W911NF- 
09-1-0248), NSF-PHY (1104546), and the Jeffress Memo- 
rial Trust (J-992). 



* Email: chuanwei.zhang@utdallas.edu 
[1] M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005) ;M. 
Dawber, K. M. Rabe, and J. F. Scott Rev. Mod. Phys. 
77, 1083 (2005); D. N. Basov et al Rev. Mod. Phys. 83, 
471 (2011); G. Catalan et al, Rev. Mod. Phys. 84, 119 
(2012). 

[2] Y. Tokura and S. Seki, Adv. Mater. 22, 1554 (2010). 
[3] T. Kimura, Annu. Rev. Mater. Res. 37, 387-413 (2007). 
[4] S.-W. Cheong and M. Mostovoy, Nature Materials 6, 13 
(2007). 

[5] R. Ramesh and N. A. Spaldin, Nature Materials 6, 21 
(2007). 

[6] W. Eerenstein, N. D. Mathur, and J. F. Scott, Nature 

442, 759 (2006). 
[7] Y. Tokura, Science 312, 1481 (2006). 
[8] G. Jin, K. Cao, G.-C. Guo, and L. He, Phys. Rev. Lett. 

108, 187205 (2012). 
[9] I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 

(2006). 

[10] H. Katsura, N. Nagoasa and A. V. Balatsky, Rev. Mod. 

Phys. 95, 057205 (2005). 
[11] M. Greiner et al, Nature (London) 415, 39 (2002). 



[12 
[13 
[14 
[15 
[16 

[17 

[18 

[19 

[20; 

[21 

[22' 
[23; 

[24' 
[25' 

[26 
[27 

[28; 
[29; 
[30; 

[31 

[32; 

[33' 
[34 

[35; 

[36; 

[37; 
[38; 

[39; 
[40; 

[41 



I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

A. B. Kuklov and B. V. Svistunov, Phys. Rev. Lett. 90, 
100401 (2003). 

L. -M. Duan, E. Dernier and M. D. Lukin, Phys. Rev. 
Lett. 91, 090402 (2003). 

E. Altman, W. Hofstetter, E. Demler and M. D. Lukin, 

New Journal of Physics 5, 113, (2003). 

S. Trotzky, P. Cheinet, S. Foiling, M. Feld, U. Schnor- 

rberger, A. M. Rey, A. Polkovnikov, E. A. Demler, M. D. 

Lukin and I. Bloch, Science 319, 295 (2008). 

J. Struck, C. Olschlager, R. Le Targat, P. Soltan-Panahi, 

A. Eckardt, M. Lewenstein, P. Windpassinger and K. 

Sengstock, Science 333, 996 (2011). 

Y.-J. Lin, K. Jimenez- Garcia, and I. B. Spielman, Nature 
471, 83 (2011). 

S. Chen, J.-Y. Zhang, S.-C. Ji, Z. Chen, L. Zhang , Z.-D. 
Du, Y. Deng, H. Zhai, J.-W. Pan. larXiv:1201.60181 
Z. Fu, P. Wang, S. Chai, L. Huang, J. Zhang, Phys. Rev. 
A 84, 043609 (2011). 

P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, 

H. Zhai, J. Zhang, arXiv: 1204. 1887 

Y. Zhang, and C. Zhang. larXiv: 1203.23891 

L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T . Yefsah, 

W. S. Bakr, M. W. Zwierlein. larXiv: 1205.34831 

M. Kohl et al, Rev. Mod. Phys. 94, 080403 (2005). 

Y. Zhang and C. Zhang, Phys. Rev. Lett. 108, 035302 

(2012). 

Y. Zhang, G. Chen, C. Zhang. larXiv: 1111.47781 

A. C. Hewson, The Kondo Problems to Heavy Fermions 

(Cambridge University Press, Cambridge, England, 

1997). 

I. E. Dzyaloshinskii, Sov. Phys. JETP 19, 960 (1964). T. 
Moriya, Phys. Rev. 120, 91 (1960). 

I. A. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 
(2006). 

C. Sen, S. Liang, and E. Dagotto, Phys. Rev. B 85, 
174418 (2012). 

S. Liang, M. Daghofer, S. Dong, C. Sen, and E. Dagotto, 
Phys. Rev. B 84, 024408 (2011). 

S. Dong, K. Yamauchi, S. Yunoki, R. Yu, S. Liang, A. 
Moreo, J.-M. Liu, S. Picozzi, and E. Dagotto, Phys. Rev. 
Lett. 103, 127201 (2009). 

M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). 

H. Katsura, N. Nagaosa, and A. Balatsky, Phys. Rev. 

Lett. 95, 057205 (2005). 

G. R. Blake, L. C. Chapon, P. G. Radaelli, S. Park, 
N. Hur, S-W. Cheong, and J. Rodriguez- Car vaj a, Phys. 
Rev. B 71, 214402 (2005). 

T. A. Corcovilos, S. K. Baur, J. M. Hitchcock, E. J. 
Mueller, and R. G. Hulet, Phys. Rev. A 81, 013415 
(2010). 

W. S. Bakr, J. I. Gillen, A. Peng, S. Foiling, and M. 
Greiner, Nature 462, 74 (2009). 

C. Weitenberg, M. Endres, J. F. Sherson, M. Cheneau, 
P. Schau/3, T. Fukuhara, I. Bloch, and S. Kuhr, Nature 
471, 319 (2011). 

J. Radic, A . Di Ciolo, K. Sun, V. Galitski, 
larXiv: 1205.21101 

W. S. Cole, S. Z. Zhang, A. Paramekanti, and N Trivedi, 

larXiv: 1205.23191 

Z. Cai, X. Zhou, C. Wu. larXiv:1205.3116l 



